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We present a method for evaluating the partition function in a varying external field. Specifically, 
we look at the case of a non-interacting, charged, massive scalar field at finite temperature with an 
associated chemical potential in the background of a delta-function potential. Whilst we present 
a general method, valid at all temperatures, we only give the result for the leading order term 
in the high temperature limit. Although the derivative expansion breaks down for inhomogeneous 
backgrounds we are able to obtain the high temperature expansion, as well as an analytic expression 

■ for the zero point energy, by way of a different approximation scheme, which we call the local Born 

, approximation (LB A). 

^ I PACS numbers: ll.lO.Wx, 05.30.Jp 

Keywords: Quantum fields, Finite Temperature Field Theory 

D 



in 



February 1, 2008 



In this Letter we discuss the evaluation of the partition function in the non-smooth backCTOund of a spherically 
fS| ', symmetric shell, a particular motivation being some recent interest in singular potentials P, ■ The delta- function 
^ profile is a useful approximation in modelling semi-transparent boundary conditions, which are naturally expected in 
many problems of physical interest. The Casimir energy has been evaluated for a scalar field in the background of 
> a spherical penetrable shell in 0|, which used the Jost function technique along with contour integral methods, e.g., 
see and the references therein. 

However, here we shall investigate the partition function for a massive, charged scalar field at finite temperature 
' with a non-zero chemical potential. We shall rely on the fact that the thermal partition function can be related (by 
, partial wave analysis) to a radial momentum integral and angular momentum sum over the phase shifts. The phase 
shift method has a simple connection with the Jost function approach . Numerous works have employed the partial 
waves technique to evaluate one loop and non-perturbative effects. Some of these are for example; the prefactor in 
bubble nucleation 0, Q , quantum effects for solitons ||| , the Casimir energy of a skyrmion and instanton effects 
^ in QCD nil Indeed, the phase shift method appears to have been first used in the context of field theory by 
Schwinger [l3| . 

^ ' Usually, the derivative expansion can be employed as an approximation to external fields that gradually vary, but 
. ^ , for the case of potentials such as a delta-function this approximation clearly breaks down. However, recent work which 
• deals with the heat kernel asymptotics of singular potentials [T3 . IT^ means that other approximation schemes can be 

J-j ' devised. We shall explicitly show how by obtaining the high temperature limit using the local Born approximation 
. 9^. , (LBA). 

In tandem we shall also apply the LBA to the calculation of the zero point energy, which gives surprisingly good 
agreement with exact results. The zero point energy must undergo regularisation; and of the many possible approaches, 
such as subtracting terms from the Born series 3i tackle the problem by subtracting the relevant heat kernel 
coefficients, enabling an analytically continued expression for the zero point energy. As discussed in |^ the intrinsically 
local heat kernel coefficients can be related to the non-local Born series. Indeed a non-local generalisation of the heat 
kernel asymptotics exists in the guise of covariant perturbation theory, which has been applied to finite temperature 
field theory in [Tij . 

As is well known, e.g., see jl^, the factored one loop effective action for a charged scalar field is — W^{ii) + 
W^{—^), where 



W'^{±fi) = ilogdet£2 



-De + - ±2ifj.-^ 
at 



(1) 
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and the parameter ^, with dimensions of length, keeps the argument of the determinant dimensionless. The connection 
between the EucUdean effective action and the canonical free energy is well known 



1^ 



(2) 



where we are ignoring any terms independent of the temperature, e.g., see We must deal with the eigenvalues of 
the operator 



A±0±(a;,r) = A„_fc(?!)±(x, r) 

where 

and the Euclidean d'Alembertian operator is 



(3) 
(4) 

(5) 



with = (Ao,0) and Aq — — i/i. It is related to the Lorcntzian □ operator by the Wick rotation —idt = dr- The 
vector Af^ is a fictitious Abelian gauge potential which represents the effect of an external charge density on the 
quantum system. 

Consider the spatial Laplacian, 



(-V^ +m^ + V{x) ) (f>{x) = 



(6) 



for a set of eigenvalues , where we shall assume a spatially dependent background coupling with spherical symmetry 
of the form 



Vir)^-S{r-R) 



(7) 



a > 0, i.e., we are considering a massive, charged scalar field in the background of a repulsive spherical shell, at some 
radius R. We have chosen the the inverse radius for our dimensionful parameter in the background delta function. 
Other choices can of course be made, see Q. If we are interested in a massless field, to = 0, a slightly different 
approach must be applied, e.g., see 0. 

After separating the eigenmodes into radial functions and spherical harmonics, the phase shift can be obtained from 
the solutions of the [d + l)-dimensional radial wave equation 



-d; - - dr + m^ + — 5{r - R) + 

r R 



u{r) = E"^ u{r) 



(8) 



with angular momentum I. The limit a — > oo implies reflecting Dirichlet boundary conditions and the partial waves 
approach should reproduce the functional determinant for the Dirichlet ball which has been studied in detail in 
(however, also see the comments in 4|). It is straightforward to show that the phase shift is (for d = 2) 



tan 5i{k) 

where the radial momentum is defined by 



^ a Ji+y2ikR)Ni+y2ikR) - 1 



k = \/ E^ — m? 



(9) 



(10) 



and Ju{x), N^(x) are Bessel and Neumann functions respectively. For such a case when the eigenvalues are not 
explicitly known it is possible to employ the phase shift method as follows. 

From a knowledge of the phase shift it is straight forward to relate it to the heat kernel by 



O /"OO 

Kit) = - / &e-('='+™')*fciV Xi'5i(fc) 
Jo , 



(11) 



e.g., see 00- We should mention that in this Letter we shall assume that there are no bound states, which is the 
case for background coupling with a > (e.g., see 0,0 for a discussion of the inclusion of any bound states). The 
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degeneracy factor = (2Z + 1) in three dimensions. There is also a free space contribution, which is encoded by the 
Cq term in the heat kernel ex pan sion. This term gives the usual constant background results, which have been well 
studied at finite temperature ^o we shall not discuss it. The heat kernel can now be used to regularise the 

one loop effective action. Defining the generalised C-function by 

C(s) = T^ /"i'"'trif(i)dt (12) 



m Jo 

The analytic continuation of ({s) then gives the one loop effective action, which is related to the zeta-function by 

w^ = w^{f,) + w''{-f,) w''{±^,)^-^c±{o) + lc±{o)loge' (i3) 

Because we are working on a static manifold, the field can be written in separable form as 4>{t,x) = e~^'^"'^(j){x), 
where w„ = 27rn//3 are the Matsubara frequencies with /3 the inverse temperature. Thus, the coincidence limit of the 
heat kernel takes the following form 

Kiit)^ E ^/ dke-^''"+'^"^'ktJ2xMk)e-^"-^'^^"' (14) 

The thermal effective action for a charged scalar field is therefore 

W^i±^i) = -- dtt'-' y - dfce-C^'+^'^'fciV x/<5i(fc)e~^""^*'^^'* (15) 

The analytic continuation can be conveniently performed by using the Jacobi-Poisson resummation formula on the 
Matsubara modes, e.g., see j2i|. After some formal manipulations we find the total thermal effective action to be 

= + W'^iti) + W'^i-^i) (16) 

In an abuse of our previous notation, see the first term is the zero point energy (/3 oo) which is ultraviolet 
divergent, hence Si{k) (see below), while the other two terms contribute to the thermal part of the effective action, 
one for each charge ±/i. Similar expressions can be found in 6, 7], however, here we have incorporated a chemical 
potential. 

For numerical purposes, the analytic continuation of the zero point energy is best performed by subtracting terms 
from the heat kernel. As t — > 0, the heat kernel in d + 1 dimensions has the asymptotic expansion 

K{t) - t-('i+i)/2 C„(r)r (17) 

where due to radial symmetry the heat kernel coefficients are only local functions of r. The leading terms, which 
cause the poles in the zeta-function, can be removed by replacing the sum over phase shifts in 1111) by (Tj 

7rCi(r)(fc2+^2)^ ^^^^^^ ^ 

xMk) = 2^ xMk) 



r(^) r(f) 

nC2{r) (fc2 + m2)^ ^^5/2(0 (fc^ + m^Y 



r(^) T{^) 



(18) 



where for distributional sources the heat kernel coefficients only contribute to surface terms [13 . ITHI |. see below. 
To obtain a finite regularised sum we only require up to C'2 for finite temperatures in three dimensions and zero 
temperature in four dimensions. However, in practice it is useful to include the C5/2 term for better numerical 
convergence or to interpolate the large k behaviour. A useful check is that the leading asymptotic behaviour follows 
that of the C5 /2 term. Indeed, as we shall show, the dominant contribution to the LBA is given by the C5 /2 term. 
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The analytic continuation of the zero point energy can be found in 0,^3 for example, and we obtain 



W°° = - 



7r/3 C. 



dk 



5/2 



2m 



+ 



Y^{21 + l)5i{k) - 2x/^Vfc2 + m2 Ci - tt C3/2 
log m i 



+ m^ 



TT C- 



5/2 



fc2 



(19) 



In Fig. 1 is a plot of the integrand (in brackets) in the above equation, given as blue solid lines for various values of 
the coupling a. The extra factor of two for the zero point energy H19|l as compared to the expression given in is 
due to the fact that we are considering a charged scalar field. Although in this Letter we do not give explicit results 
that depend on apart from l|16|) . the more general expression is required if one wishes to comment on things such 
as Bose-Einstein condensation, which we hope to report on in the near future. 

For the spatial wave equation © with a potential of the form defined in lO the heat kernel coefficients are [T^ [Tsj : 



Ci{r) 



C2{r) 



(47r)('^+i)/2 



(47r)('^+i)/2 
lo^ 



a 
R 



(4^)(d+i)/_ 



R 



^5/2(0 - (4^)(d+l)/2 



32 i?4 



i?2 



where the volume of a spherical shell in c? + 1 dimensions is given by 



2^(d+l)/2^d 



r(^) 



= AttR-" 



d=2 



(20) 
(21) 

(22) 



In fact, the leading order contribution comes from the C5/2 term in H19|l . up to the renormalisation scale dependence, 
£. As can be seen in Fig. 1, for example, the choice a = for various values of x gives surprisingly good agreement 
with the integrand of ((T^ . That is, the the zero point energy is dominated by the term 



dk 



TT a 



5/2 



fc2. 



Tip C. 



5/2 



2m 



(23) 



for which the integrand of the left hand side (in brackets) we have plotted in Fig. 1 as red dashed lines. Clearly, 
the accuracy of this approximation improves as the area of the integrand in (|19|l becomes larger and the oscillatory 
behaviour for the large k part has less effect. Including higher order corrections, such as the C3 term, should improve 
the accuracy for smaller values of the area of the integrand in H19|l . 
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z = kR 



FIG. 1: A plot of the integrand in p9|l . blue lines, and the contribution from the integrand in H'Z'A^ . red dashed lines, for values 
of increasing a = 1.5, 2.25, & 2.75, in units of (mR)"^. For clarity the z axis has been shifted below the origin by a small 
amount. 

For large temperatures numerical evaluation of W^{±iJ,) in H16|l becomes difficult and it is useful to have analytic 
expressions valid in the high temperature limit, f3 0. We can achieve this by using the LBA, as we have already 
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done for the zero point energy. The high temperature expansion can be obtained in various ways, however, we shall 
simply adapt the method used in [l3|. First, let's rewrite the thermal effective action explicitly in terms of k as 



(24) 



„/3(fc2+m2)l/2 ^^p^ _ ^ 



Our approach will be to substitute the sum over phase shifts 

n (k 



2^ xA{k) ^ ^ ^td+i-2n+2^ (25) 



into our expression for (|^ . i.e., W^{±r) — W^{±r) 

Wni^r)^ p^ d+l-2»+2 y^ ^ ^^g(a;2+m2)i/2 g±™ _ ;L (^6) 

where we have made the change of variables x = (3k, fh = fim and r = /z/m. 

In this Letter we shall only consider the leading order contribution to the high temperature expansion (for c? = 2), 
which corresponds to n = 1, i.e., 



Wf(±7-) = -/3-^— i- / rrdx— — -— — — (27) 



The small m expansion of the above integral is well known 18] and leads to the result 



<(±r) = ^52(m,±r) (28) 
where the function 52 ("^, J') is defined in Appendix A of |l8j |. Then, it is simple to show that 



g2{'m, r) + 32 (™, — f) — -z 2to + mr In 



(1 + r) 



m2(l — r^) 



©(m^) (29) 



.(1-0 

and thus, the leading order contribution for small ffi to the thermal part of the effective action is 

where the second equality is specific to a spherical shell. Thus, for the free space field theory, i.e., in a constant 
background, the leading order contribution to the free energy is oc e.g., see whereas the effect of a 

repulsive delta- function potential is to contribute to leading order oc (3~^. 

The method we have discussed is sufficiently general such that we can choose any spherically symmetric background 
potential, smooth or inhomogeneous, the only requirement being a knowledge of the heat kernel coefficients. Of 
course if they are not known then the heat kernel coefficients can be derived by taking the local part of the Born 
approximation, which is an expansion in powers of small coupling of the phase shift. This is also true for the case 
when there is no analytic expression for the phase shift, where recourse can be made to the WKB method P. IT^. 

Furthermore, the finite temperature expression (|24|l is a general expression valid at any temperature, which requires 
numerics for the general case. However, it is instructive to have analytic expressions for a charged scalar field in the 
high temperature limit, for example to consider Bose-Einstein condensation. In this Letter we showed how to obtain 
the leading order term. Indeed, next to leading order corrections will depend on the chemical potential, fi. The low 
temperature expansion can also be obtained in a similar way. 

The local Born approximation (LBA) was also applied to give an approximate expression for the non-thermal part 
of the one loop effective action, i.e., the zero point energy. We employed zeta- function and heat kernel methods to 
subtract the divergences inherent in the zero point energy, this being one of the many commonly used subtraction 
procedures. 

As well as for delta-function type potentials one could consider step function profiles, which are useful in modelling 
many physically interesting situations, such as instantons in bubble nucleation For such a case it would also be 
possible for the chemical potential, /x, to vary with the same radial profile as the mass, e.g., a step- function with 
profiles ^^9{r — R) and rn?9{r — R). However, this is a considerably more complicated set up, which we shall report 
on in forthcoming work, as well as other related issues. 
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